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ABSTRACT 



High precision radial velocity (RV) measurements in the near infrared are on high demand, especially in the context of exoplanet 
search campaigns shifting their interest to late type stars in order to detect planets with ever lower mass or targeting embedded 
pre-main-sequence objects. 

ESO is offering a new spectrograph at the VLT - CRIRES - designed for high resolution near-infrared spectroscopy with a comparably 
broad wavelength coverage and the possibility to use gas-cells to provide a stable RV zero-point. 

We investigate here the intrinsic short-term RV stability of CRIRES, both with gas-cell calibration data and on-sky measurements 
using the absorption lines of the Earth's atmosphere imprinted in the source spectrum as a local RV rest frame. Moreover, we also 
investigate for the first time the intrinsic stability of telluric lines at 4100 nm for features originating in the lower troposphere. 
Our analysis of nearly 5 hours of consecutive observations of MS Vel, a M2II bright giant centred at two SiO first overtone band-heads 
at 4100 nm, demonstrates that the intrinsic short-term stability of CRIRES is very high, showing only a slow and fully compensateable 
drift of up to 60m/s after 4.5 hours. The radial velocity of the telluric lines is constant down to a level of approx. ±10m/s (or 7/1000 
of one pixel). Utilising the same telluric lines as a rest frame for our radial velocity measurements of the science target, we obtain a 
constant RV with a precision of approx. ±20 m/s for MS Vel as expected for a M-giant. 

Key words. Instrumentation: adaptive optics - Instrumentation: spectrographs - Methods: observational - Techniques: radial veloc- 
ities - Techniques: spectroscopic - Stars: AGB and post-AGB 



1. Introduction 

CRIRESE] is a newly commissioned adaptive optics (AO) - fed 
spectrograph at UT1 (Antu) of the Very Large Telescope on 
Paranal, Chile. The spectrograph provides a resolution of up 
to 100000 at 960-5200 nm with an instantaneous wavelength 

coverage of ~ /1/ 50 on a 4096x512 pixel detector mosaic 

dKaufl et al.L 12004 120061 ). A MACAC0 adaptive optics system 
dArsenault et al.L 120041) is providing diffraction limited images 
at the entrance slit of CRIRES with Strehl ratios of up to 65% in 
the K band, thus roughly doubling the throughput under reason- 
able seeing conditions for the nominal 0.2" slit whe n compared 
to seeing limited performance (Paufiqu e et aUl2006l) . 

The spectrum, being either a full order at the shortest wave- 
length in the J band and about 1/5 of an order at the longest 
wavelength in M band, is imaged in long slit mode onto a mo- 
saic of four Aladdin III detectors with a gap of ~ 250 pix be- 
tween the chips. Thus, the instantaneous wavelength coverage 
is much shorter than for optical cross-dispersed spectrographs. 
However, the gain in efficiency for RV measurements, primar- 
ily dependent on spectral coverage and S/N, will shift in favour 
of CRIRES as soon as the effective SED of the science target 
peaks longwards of 900 nm where other optical spectrographs 



1 CRIRES stands for CRyogenic InfraRed Echelle Spectrograph 

2 MACAO stands for Multiple Application Curvature Adaptive 
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commonly used for RV measurements, are either not operating 
(e.g. HARPS) or loose dramatically in sensitivity (e.g. UVES). 

For precise RV measurements the spectrograph has to be 
well characterised since we deal with various sources of un- 
certainties that could induce wavelength shifts in the measured 
spectra. Since CRIRES is a cryogenic, thus actively stabilised 
spectrograph, temperature drifts are small but certainly not negli- 
gible. The same is true for vibrations and flexure, especially due 
to changing thermal gradients in the structure. Being mounted on 
a Nasmyth platform of the VLT, the system is intrinsically very 
stable but vibrations induced by the cryogenic cooling system 
can be suspected for causing instability. 

A second effect is the imperfect repeatability of dispersive 
parts on moving functions suffering slip-stick effects: the prism 
(pre-disperser), the intermediate slit and the grating. It is forseen 
that functions are left untouched when the same instrumental 
setup is used in subsequent observations. However, for obser- 
vations that span days to months the repeatability of the setup 
could be an issue, if it can't be well calibrated. 

Last but not least a long-slit spectrograph suffers the natural 
problem of alignment uncertainties of the source with respect 
to the slit centre. For spectrographs fed by an adaptive optics 
system the situation is worse in some aspects: 

1. The FWHM of the source is usually smaller than the slit 
width. 

2. The alignment of the adaptive optics in respect to the en- 
trance slit of the spectrograph suffers instabilities, requiring 
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an active slit-viewer guiding in order to keep the source cen- 
tred on the slit. 

3. The algorithm for slit-viewer guiding can only work on the 
wings of the stellar PSF, since its core is masked by the slit. 

4. The refraction caused by the Earth's atmosphere has to be 
measured and compensated simultaneously as observations 
are performed at three different wavelengths: in the optical 
for the AO wave front sensor, in J, H or K band for the slit- 
viewer (used for active guiding) and at the respective wave- 
length of the observation, in our case at 4100 nm. 

All these effects will cause a degraded stability and a possible 
drift of the wavelength scale with time. We thus aimed to charac- 
terise the system and the means to calibrate it and test the on-sky 
performance in terms of RV stability during a single night. 

For CRIRES a RV precision of 15 m/s would require point- 
ing and tracking with a precision of one milli arcsec (mas), given 
that one pixel in dispersion direction (1.5 km/s) equals 100 mas 
at the entrance slit of the spectrograph. The tracking precision is 
limited by the amount of light in the wings of the PSF on which 
the guiding algorithm of the slitviewer operates. Tracking the 
core of the PSF on the seeing-limited wings is a general limi- 
tation. Even for bright sources, the performance of the guiding 
algorithm is not sufficient to stabilise the source in the slit with 
a precision better than approx. 25 mas (approx. 200 m/s), which 
makes a reference source (telluric or gascell absorption lines) in- 
dispensable when aiming for high precision RV measurements. 

2. Observations 

For our on-sky test we were looking for a bright target which 
should present no significant R V signal with in a few hours. 
We chose the M2II star MS Vel (lEgrett 1 19801). a rather brigh t 
(Ks = 3.7 mag (DENIS), F^^f 16Jy (lEgan et all 120011 
semi-regular variab le star dAdelmanl 200 lb with a proposed pe- 
riod of about 360d (Ma karov et al.L 119941) . We have re-analysed 
HIPPARCOS and TYCHO photometry simultaneously cover- 
ing about 3 years and about 5.5 years w orth of data from the 
ASAS3 V band survey dPoi manskil |2002|) using a Lomb-Scargle 
Periodogram to identify variability in MS Vel. All three data sets 
show a strong period of P = 360 + 5d with a somewhat larger 
value (P ~ 363.5 d) for the HIPPARCOS and TYCHO epoch of 
1990 - 1993 and a shorter value (P ~ 355 d) in the full epoch 
covered by the ASAS data (2000 - 2006). This dominant period 
with an amplitude of about AV ^ 0.8 mag is accompanied by 
shorter periods of 18, 26 and 28 days with much smaller ampli- 
tudes in the earlier data which disappear however in the much 
better sampled data from the later epoch. Conversely, periods of 
155 and 185 days, present in the ASAS data can not be found 
in the HIPPARCOS and TYCHO data. All these periods have 
false-alarm-probabilities of much less than 0.1% and are there- 
fore real, even though not stable. We see that MS Vel is variable 
on various time-scales of weeks to a about a year. However, it 
is very unlikely that the object shows a significant RV signal in 
only a fe w hours, given its large ra dius. As predicted by 3D sim- 
ulations dFrevtag & Hofnen, 12008), typical photometric periods 
for early to mid M-type g iants are at the ord er of a several dozend 
to a few hundred days ( Perc y et all 120081) . Radial velocity sur- 
veys of K-type giants also revealed periods at the same timescale 
dHekker et all 12008b . Short-term oscillations on M-type giants 
are rarely observed an d still show periods in excess of one day 
dKoen & Lanevil2000l) . 

We have observed MS Vel in the night of 2007, Feb 02 from 
03:23 to 07:56 UT. Observations were obtained with the adap- 



tive optics loop closed and locked on the science target. Strehl 
ratios of more than 50% could be achieved during the majority 
of the program, delivering most of the encircled energy in the 
diffraction limited core with a FWHM of approx. 100 mas. The 
DIMM seeing was 1" on average, with very good conditions 
(0.65"-1.0") for the first three hours and slightly worse condi- 
tions (1.0"— 1.5" ) for the last 90 min. A slit width of 0.4" was 
chosen, providing a nominal resolution of R^ 50 000. The total 
flux in the spectrum shows a weak dependence on the seeing, 
with a degradation of about 20-30% in flux at the end of the 
campaign, as expected when the atmospheric conditions wors- 
ens and the Strehl ratio drops, matching the predictions of the 
exposure time calculator (ETC) for our case very nicely. 

We chose a setting at A^f = 4029.9 nm to cover two 
SiO band-heads at 4043 nm ( 28 Si 15 0, v=3^1) and 4084 
nm ( 28 Si ls O, v=4— >2). The range of the SiO-overtone at 
4000^4100nm is a well understood but still rich line system 
dAringer et al.L 1 19991) in a region where photospheres are other- 
wise generally clean while the Earth's atmosphere provides for 
a rich line system of N2O, originating mainly from the lower 
troposphere. 

An AAA-BBB-BBB-AAA nodding pattern with individual 
integration times of 10 s per frame was used. The nodding am- 
plitude was set to 10". The average S/N in each single spectrum 
was found to be approx. 200. 

After each cycle the instrument switched to calibration mode 
to take a single, short spectrum of the N2O gas-cell (10 hPa, 
140 mm length), illuminated by a blackbody involving an 
Ulbricht sphere. The slit width was held constant for the science 
and calibration measurements. Since the light from the gas-cell 
is a spatially flat source, the entrance slit of CRIRES is evenly 
illuminated and the 0.4" slit delivers a resolution of approx. 
50000. No change in illumination is expected throughout the 
length of the program, making these calibration measurements a 
suitable reference for the intrinsic stability of the spectrograph, 
independent from slit centring issues of the actual science target. 
After completing such a set of observations we re-acquired the 
object in the slit-viewer, closed the AO loop again and started a 
new cycle of on-sky observations. This procedure was repeated 
25 times. 

During this series of observations, the dispersive elements 
of CRIRES - the prism (pre-disperser), intermediate slit and the 
grating - where electronically de-mobilised and held in fixed po- 
sitions. This precaution was taken to avoid accidental reposition- 
ing of functions. Thus, the observations of the calibration source 
were only affected by temperature drifts and flexure. Besides 
these non standard setup the observations made use of ESO ob- 
servatory standard procedures. 

Both, the signal-to-noise ratio and the spectral resolution 
are of critical importance for RV studies. While both quanti- 
ties are closely connected to the slit width of regular long-slit 
spectrographs, the situation for CRIRES is more complex. Since 
CRIRES works with an underfilled entrance slit, the achievable 
resolution is not strictly depending on the chosen slit width but 
is mainly governed by the FWHM of the stellar PSF in the slit, 
hence, also depends on the AO performance and the wavelength 
of the observation. The troughput is only marginally dependant 
on the slit width in the regime of potentially high Strehl ratios 
in K, L, and M band, reaching over 80% at the longest wave- 
length under favour able conditions (see the CRIRES ETC and 
iPaufique et alll2006l) . The typical diffraction limited FWHM of 
the VLT at 4000 nm (~0.1" ) is much smaller than the minimum 
slit width and the spectral resolution is limited by the optics of 
the spectrograph and the stability of the source in the slit during 
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an integration rather than by the width of the profile alone or the 
width of the entrance slit. The true resolution achieved with the 
0.4" slit is thus likely closer to 100 000 than to its lower limit of 
50000, defined by the slit width. Closing the slit to 0.2" would 
only marginally improve the resolution by basically masking 
some light from the seeing limited halo of the PSF. Likewise, 
the troughput of the 0.4" slit is only about 20% higher than for 
0.2" slit in our case. 

However, choosing a wide slit further decreases the amount 
of reflected light available for the guiding algorithm of the 
slitviewer and further increases the chance of misplacements and 
instabilities of the source position in the slit. Since it is expected 
that future RV programs conducted with CRIRES will also ob- 
serve rather faint targets under low Strehl ratios, a regime where 
a large slit width will more notably improve the efficiency, we 
nevertheless decided to use the 0.4" slit to investigate these ef- 
fects also in this first study. 

3. Data Reduction 

3.1. Basic Treatment and Spectral Extraction 

For data reduction the raw data were divided by a flat-field 
frame, mainly to correct for a fixed pattern effect of the detec- 
tor array and to normalise the continuum to a satisfactory level. 
To subtract the background emission we subtracted the mean of 
the three frames from the next nodding position from the current 
frame, leaving a flat sky background around the source spectrum. 
The spectra were than extracted using the optimum extraction al- 
gorithm of the CRIRES data reduction pipeline. 

Spatially extended sources, such as the sky background as 
well as the spectra of the gas-cell in the respective calibration 
frames have been extracted after flat-fielding of the raw frames 
and integrating over a spatial range of 100 pixels between the 
two nodding positions of the source. We thus sample the same 
region that is covered by the source spectrum itself, without in- 
troducing offsets caused by slit curvature, detector residual tilt, 
and distortion, as we will discuss in Sect. [4] 

After this treatment we have 300 individual spectra of the 
object and of the sky background emission at our disposal. In 
addition we have 25 spectra, one for each cycle of observations, 
of the N2O gas-cell. Hence, we have four systems of spectral 
lines that we can analyse here: 

1. the photospheric lines of the target, MS Vel, foremost the 
two SiO band-heads and various unidentified lines, 

2. the telluric lines in absorption, imprinted in the spectra of the 
target, 

3. the telluric lines in emission as a spatially extended back- 
ground source along the whole slit (note: due to Kirchhoff 's 
law the atmospheric absorption lines become inverted into 
emission lines when viewed against the night sky) and fi- 
nally 

4. the N2O lines from the gas-cell, also as a spatially extended 
source since the instruments flat-field source was used to il- 
luminate the gas-cell. 

All four line systems can now be analysed by means of cross- 
correlation to measure the zero-point shift of the wavelength in 
each of the spectra. 

3.2. Methodology 

In case the spectrum of the star is taken through a gas-cell, 
the classical approach to measure radial velocities is to cross- 



correlate the product of two template spectra (normally the stel- 
lar spectrum and the spectrum of the gas-cell) with the actual 
measurement. This requires template spectra of high quality. 
Unfortunately no template spectra for a M2II giant at 4100nm 
of sufficient quality was available for us. A test with a synthetic 
spect rum from the library of PHOENIX ([Hauschildt & Baron, 
2005, and references therein) spectra^, as displayed in Fig. Q] 
did not lead to consistent results when used as a template for the 
cross-correlation analysis, since the line lists used to compute the 
synthetic spectra are not precise enough and show notable shifts 
of the different species present in late type giants at 4100nm. 

For the telluric lines however, we have access to quite pre- 
cise synthetic spectra of the Earth's atmosphere. Using a stan- 
dard atmospheric model for the Paranal site and adopting a 
typical humidity for the time of the observation we can con- 
struct a synthetic spectr um at any given wa velength using the 
FASCODE algorithm dClough et all Il98ll |1992|) . a line-by- 
line radiat ive transfer model for the Earth's atmosphere, and 
HITRAN dRothman et all 120051) as a database for molecular 
transitions. A spectrum of the Earth's atmosphere in our wave- 
length window is shown in Fig. Q] for the telluric lines in ab- 
sorption as seen against the stellar continuum. In Fig. [2] we show 
the same lines in emission (line radiation against dark sky back- 
ground). It is worth noticing that the predominant species in the 
observed spectral range is in fact N2O, the same species as in our 
gas-cell. N2O shows nearly regular spaced lines of the f-branch 
on chips two, three, and four and of the /^-branch on chip one. 
Additional minor contributions arise from some isolated water 
vapour lines. These lines form in lower, thus warmer parts of 
the atmosphere and are stronger in emission than in absorption 
when compared to the neighboring N2O lines. They are thus 
more prominent in the sky emission spectrum (Fig. [2]) than in its 
'inverse' counterpart, the telluric absorption spectrum (Fig.[TJ. 

A more detailed discussion of the properties of the telluric 
lines and the performance of using them as wavelength standards 
will be given in a later paper. 

Similarly to the telluric features, we can compute the spec- 
trum of our gas-cell, a N2O spectrum at room temperature and 
~10 hPa pressure. The result is shown in Fig. [3] These transitions 
are calibrated in f requency with heterody ne techniques against 
the time standard dMaki and Wellsl |1992|) . so that they provide 
for an absolute precision, traceabl e to the t i me sta ndard, not yet 
available for optical spectroscopy dLi et all l2008). 

3.3. Wavelength Calibration 

For a RV measurement we need to determine the wavelength 
of each pixel to a high precision. The gas-cell delivers a N2O 
absorption spectrum of high S/N with nearly evenly spaced lines 
being distributed in high density on all four chips (see Fig. [3). 
Using the fore-mentioned synthetic spectrum of our gas-cell, we 
fitted a quadratic solution independently to each of the four chips 
to match the measured spectrum with the synthetic spectrum, 
thus providing the zero point and dispersion for the four chips in 
the first recorded N2O spectrum during our measurements. 

3.4. Cross-Correlation 

Given the lack of a proper template spectrum for the science 
target, we use the following approach to determine the radial ve- 
locity of each of the line systems as outlined above: (1) For the 
photospheric lines of MS Vel we mask the telluric lines for trans- 

3 kindly provided by Peter Hauschildt, Hamburg 
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missions below 95% and cross-correlate each of the 300 spectra 
against the first one obtained. (2) The wavelength shifts of all 
other line systems - the telluric absorption lines in the spectrum 
of MS Vel, the sky-background emission lines and the N2O ab- 
sorption lines - are computed by cross-correlating our measured 
spectra against the respective synthetic model spectrum. 

We also tested the methodology used for the photospheric 
lines of MS Vel (cross-correlation of all spectra against the first 
one) for the N2O lines and the telluric lines. We obtained the 
same results as for the cross-correlation with the synthetic tem- 
plates, but find an enhanced scatter due to the lower signal-to- 
noise when cross-correlating to a measured spectrum instead of 
against a noise-free model. 

4. Results and Discussion 

The 25 calibration spectra of the N2O gas-cell are distributed in 
time over the full length of our experiment. These spectra probe 
the intrinsic stability of the spectrograph, since the illumination 
of the gas-cell (and thus of the entrance slit) is spatially flat and 
we do not trace any slit alignment effects. Moreover, the tem- 
perature of the gas-cell, even though not actively stabilised, is 
certainly not variable to a level that could induce line broaden- 
ing and line-shifts in the N2O spectrurrQ. Hence, we adopt this 
spectrum as a true rest frame. As can be seen in Fig. [4] the spec- 
tral lines of the gas-cell show a slow drift in wavelength that 
culminates in a RV offset of ~ 60m/s after 4.5 hours. The drift 
is smooth but non-monotonous and is very similar for all four 
chips leading only to a minimal degree of dispersion change, as 
can be seen by the slightly diverging RV trends of the four chips. 
Such a change in dispersion can be assigned to a small drift in 
the temperature of the grating. At the time of our measurements, 
the grating was only passively stabilised and the temperature 
sensor at the grating had a limited resolution of only ~ 0.1 K. 
Moreover, the sensor was not part of an active temperature sta- 
bilisation of the grating. Following the basic grating equation 
for echelle spectrographs, the angular dispersion changes lin- 
early with the groove spacing of the grating (induced e.g. by 
thermal expansion of the grating material) when the setup is oth- 
erwise held constant. Adopting a typical expansion coefficient of 
dl/l - 4 ■ 10 6 Kr 1 for the grating material, a temperature drift 
of only 50 mK over 4.5 hours could have caused the observed 
drift of 60m/s. Such a drift would have easily passed the sensor 
unnoticed. This shortcoming has been fixed by now. The grat- 
ing is now actively stabilised to 1 mK. The expected wavelength 
stability is thus approx. 1 m/s. 

This drift would be seen in any spectral line that is intrinsi- 
cally stable and used as a local rest frame and can thus be com- 
pensated for. Ideally the gas-cell should be used simultaneously 
while measuring the spectra of the science target, the absorp- 
tion lines of the gas being imprinted into the source spectrum. 
However, the chosen gas has to match the spectrum of the object 
and the gas should not be present in the Earth's atmosphere to 
avoid any overlap or mixing of the spectrum of the Earth's at- 
mosphere and the spectrum of the gas-cell. At the time of our 
experiment we had only one gas-cell ready and its content, N2O, 

4 During our observations the external temperature as logged by 
the observatory dropped from 14.8 to 13.3 C. The gas-cell has some 
thermal inertia, but will ultimately follow the external temperature. 
Assuming a ty pical pressure i nduce d frequency shift of 100-150 MHz 
for 1000 hPa jDeming et all 119871) . even the maximum temperature 
change for Paranal (-5 < T am b < 20 C) will not result in a shift of 
more than 150 kHz of the line frequencies, equivalent to a velocity er- 
ror of 50 cm/s. 



is also highly abundant in our atmosphere. Hence, we decided to 
use the gas-cell only in calibration mode, illuminated by the flat- 
field source, to trace the internal stability of the spectrograph 
and use the atmospheric lines, seen in absorption in the target 
spectrum, as our local rest frame. 

Conversely, we can use the atmospheric lines in emission to 
trace the behaviour of these lines - to check for their stability. 
Here the same argument holds as for the N2O lines: The atmo- 
spheric lines in emission are a spatially flat source unaffected by 
any slit misalignments. The frequency of the atmospheric lines 
can be systematically shifted if there were a net radial veloc- 
ity introduced by wind or turbulence. Moreover, these lines will 
be affected by the slow drift of the spectrograph itself, measured 
with the N2O gas-cell lines. In Fig.[5]we show the RV trend of the 
atmospheric lines in emission. The measured RVs from 12 spec- 
tra taken within one nodding cycle are averaged to a single data 
point. The standard deviation of the mean is shown as the formal 
lcr error of each data point. We over-plot the drift of the spectro- 
graph and subtract this drift from the RV points. Combining the 
results from all three chips where atmospheric lines are present 
results in a flat RV distribution over time with residuals (r.m.s.) 
in the order of < 10 m/s which is also well within the formal 
uncertainties derived for each data point. Thus, we conclude at 
this point that the atmospheric lines are well suited as a refer- 
ence to trace the RV of the science target. We also note here that 
the spectral resolution of the emission lines is limited by the slit 
width to R~50 000 while the resolution of the object spectrum is 
more dependent on the true FWHM of the PSF in the slit, thus 
providing a higher resolution. 

Finally we measure the RV of the photospheric lines of 
MS Vel and in parallel the atmospheric lines in absorption, as 
outlined above, that is by cross-correlating the masked spectra 
of the source against the first measured spectrum and the un- 
masked spectra against a template of the atmospheric absorption 
lines. Both line-systems show a strong near-linear trend with 
time, peaking at ~ 1.5km/s after 4.5 hours. This trend is most 
likely caused by residuals from the correction of the chromatic 
refraction (see point 4 in Sect. Q]) In addition, the nodding pat- 
tern induces sudden jumps of up to 500 m/s (1/3 pixel) between 
the two nodding positions due to the slit curvature and a slight 
rotation of chips one and two in the focal plane. 

Subtracting the RV of the atmospheric absorption lines (that 
suffer the same displacements) from the RV measured in the 
masked source spectrum cancels out most of the effects of the 
nodding displacement and the source misalignment due the re- 
fraction correction. A linear trend is left, in the counter direc- 
tion of the trend observed in the individual spectra of the source, 
peaking at 300 m/s after 4.5 hours. The results at this point of 
the analysis are displayed in Fig. [6] Again, we combine here the 
12 spectra from each of the 25 nodding cycles to one single data 
point. 

The remaining trend can be modeled by the barycentric mo- 
tion and rotation of the Earth. Subtracting this trend leaves resid- 
uals in the order of approx. + 20 m/s. The residuals appear not 
noise-like, even though the formal uncertainties of the RV points 
would account for most of the amplitude of these residuals, but 
show a rather smooth and semi-periodic behaviour. These resid- 
uals are not correlated with the residuals in the atmospheric 
emission lines (see Fig. [5j making it unlikely that the variabil- 
ity of atmospheric lines , e.g. by wind drifts, are responsible for 
the residuals in the data. It is also very unlikely that the source 
itself shows a RV variability within one or two hours, given the 
large radius of the giant star. Hence, we assign the residuals to 
a general noise floor limiting the RV precision to this level, be- 
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ing most likely due to the data analysis rather than intrinsic to 
instrumental instabilities. 

We note, that our analysis does not take changes in the line- 
profile into account. This is an essential part of the analysis when 
aiming for very high ra dial velocity precision dValenti et al.L 
119951; iButler et all \l 996) using in-situ measurements through a 
gascell. The method involves a template spectrum of the gas- 
cell with much higher resolution than the one used at the actual 
science measurements. This template spectrum is used to decon- 
volve the measured lines to determine the instrumental profile. 
We did not dispose of such a spectrum of our N2O gascell, nor 
do we have such a template for the telluric lines. Hence, this 
shortcoming in our analysis can easily account for the remain- 
ing residuals. Moreover, the masking technique to separate the 
telluric and stellar lines in the individual spectra is not perfect 
for highly blended line systems, as is the case here. 



5. Conclusions 

We have performed the first radial velocity stability test of 
CRIRES with 4.5 h of on-sky measurements of the giant MS Vel 
at 4100nm. The dispersive parts of the spectrograph have been 
electronically stabilised, thus decoupling the problems of func- 
tional repeatability, spectrograph stability and tracking errors. 
Besides this non-standard settings we used the spectrograph as 
offered to the community, which demonstrates that even at this 
very early stage of its life, CRIRES is a stable instrument, ready 
to measure radial velocities with high precision including a po- 
tential for even higher performance when using an appropriate 
gas as local RV zero point. 

We could demonstrate that CRIRES shows only a slow drift 
in its intrinsic radial velocity zero point of ~ 60 m/s in 4.5 h. This 
drift is fully compensatable using a N2O gas-cell in calibration 
mode or telluric lines as a local rest frame. Moreover, if our as- 
sumption holds that the drift is due to a small change in temper- 
ature of the grating (by approx. 1/20 K), it will not appear in any 
later measurements as the grating is now actively stabilised to 
approx. 1 mK. 

The telluric lines in emission, observed with a resolution of 
R~50 000 are stable in RV to a level of < + 10 m/s and are there- 
fore suitable as a local rest frame for high precision RV mea- 
surements. For the relative radial velocity of our science target, 
MS Vel, we could recover the 300 m/s drift induced by rotational 
and barycentric Earth's motion. The remaining residuals (rms) in 
the stellar spectra are < + 22 m/s. 

We can conclude that the telluric lines imprinted in the sci- 
ence object are not a limiting factor for precise RV measure- 
ments in the near infrared but can instead be used as a substitute 
for a gas-cell providing the necessary rest frame in cases where 
no suitable gas can be found or the interesting stellar features 
are in regions of high telluric line density. First sc i entific results 
using this method are presented in Huelam o et alj d2008l) . 

The influence of mismatches between the diffraction limited 
FWHM of the PSF and the typical slit width in combination with 
dynamic (i.e. tracking induced) misalignments of the PSF in re- 
spect to the slit centre are substantial. However, these problems 
can be compensated when telluric or gas-cell lines are measured 
simultaneously with the source and used as a local rest frame. 
This method is well known and widely used in the optical, where 
Iodine lines are used as the local restframe. Given the high rate of 
close blends of telluric lines and photospheric lines from our sci- 
ence target and taken the limitation of our analysis into account, 
we see a clear margin for a gain in performance. Especially the 



use of a gas-cell, well characterised by a high resolution labo- 
ratory spectrum (R» 100 000), should provide a better local rest 
frame when aiming at RV precision below «20m/s. The gas has 
to exhibit a high and stable line density in regions where the sci- 
ence target shows also a reach spectrum. To avoid blends with 
telluric lines, rare isotope enhanced gases could be used which 
would elegantly avoid problems with blending telluric features. 

In the example presented here - as for all conceivable RV 
applications of CRIRES - rotational vibrational transitions of 
molecules in the vibrational ground-state are used. These gases 
behave very close to ideal gases. This implies that the density of 
molecules in the beam is invariable and relatively easy control- 
lable. This situation is fundamentally different in the Iodine gas- 
cells which are being used for radial velocity work. There iodine 
vapour is used in a heated cell and condensation on cooler parts 
can result in strong changes of the vapour pressure. Therefore, 
in our case it is safe to assume, that the transition frequencies 
are very close to those, which have been determined by NIST 
using heterodyne-techniques relative to the Caesium time stan- 
dard. The spectra are relatively simple and lines which may con- 
sist of a blend can be positively excluded. The only conceivable 
process which may change transition frequencies in the gas-cell 
in the case of an ideal gas is pressure shift. The pressure may 
change due to the external temperature, which will be monitored. 
It has been estimated that this effect is less than lm/s. A more 
serious problem is pressure shift b y air, in case the gascell has 
a leak. Fig. 2b in lGlenar et alj d!988l) gives an example. A PH3- 
gascell line in the presence of a leak is measured relative to a 
gain-stabilised C02-laser in a heterodyne setup. The kinks in 
the PH3-trace correspond to refillings of the cell. To control this 
problem, the pressure in the CRIRES gas-cells will be monitored 
as part of regular operations. 

A gascell filled at low pressure with species, for which NIST 
has done a heterodyne reference measurement relative to the 
time standard (Caesium-clock), in principle provides for an ab- 
solute secondary frequency standard. To that end CRIRES mea- 
surements are linked in an absolute and traceable way to the time 
standard. This is, strictly speaking, not an actual experimental 
proof, that this indeed allows for a "cm/s"-type calibration of 
the spectrograph. Indeed e.g. a spurious blend of gas-cell lines 
with weak telluric lines may produce artificial shifts exceeding 
this value. Experts in the field indeed stress that the frequency 
calibration of spectrographs breaks down at a precision of 10 10 
and only heterodyne techniques can provide for this precision 
(Theodor Haensch, priv. communication). 

We finally note that this paper aimed at demonstrating the 
short-term stability of CRIRES and the implications when ob- 
serving with an underfilled entrance slit. The long-term stability 
of the instrument, especially important for RV studies, will be 
considered in a later paper. 
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Fig. 1. Spectrum of MS Vel and the Earth's atmosphere on the four chips of CRIRES, displayed in four viewgraphs. Shown in each 
viewgraph: measured spectrum (top, black line), telluric model (center, red line,shifted by -0.4) and PHOENIX model spectrum for 
MS Vel: T e //=3500K and logg=0.5 (bottom, green line, shifted by -0.8). See text for details. 
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Fig. 2. Spectrum of the Earth's atmosphere in emission (sky background) on the four chips of CRIRES, displayed in four viewgraphs. 
The relative strength of the emission bands follows the behavior of the N2O absorption bands shown in Fig. [3] The relative intensities 
on chips 2, 3, and 4 are scaled up by factors of 7%, 15% and 150%, respectively, relative to chip 1. Over-plotted is a synthetic 
emission spectrum (red line). Note that the isolated H2O lines appear much stronger in emission than they appear in absoption in 
Fig-Q] See text for further details. 
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Fig. 4. Apparent radial velocity of the N2O gas-cell line system (see Fig. [3} displayed separately for each chip) observed in cal- 
ibaration mode (the gas-cell being illuminated by the flat field lamp). Note the change in dispersion with time which manifests as 
divergent trends in the RV measured on the four chips. See text for details. 
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Fig. 5. Radial velocity of the telluric emission lines (see Fig. 13 and the N2O gas-cell line system (see Fig. [3]) showing the instru- 
mental drift during the measurement. In the subpanels the difference between both measurements is shown in terms of observed 
- computed (O-C), first separately for each of the three chips analyzed here, finally for the mean of all three chips. The standard 
deviation of the mean was used to derive formal error bars. Residuals are the final uncertainty for the RV stability of the telluric 
lines. See text for details. 
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Fig. 6. Radial velocity of the science target - MS Vel - after subtraction of the radial velocity of telluric absorption lines imprinted 
in the same spectrum (see Fig. Q}. The linear trend is modeled by the barycentric and rotational velocity of the Earth during the 
observation. In the subpanels the data after subtraction of the model is shown in terms of observed - computed (O-C), first separately 
for each of the three chips analyzed here, finally for the mean of all three chips. The standard deviation of the mean was used to 
derive formal error bars. See text for details. 



